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Abstract 

Calculations of the ^'^He + ^^Si total reaction cross sections at intermediate energies 
O . are performed on the basis of the Glauber-Sitenko microscopic optical-limit model. The 

I target-nucleus density distribution is taken from the electron-nucleus scattering data, 

■ and the ®'^He densities are used as they are derived in different models. The results 

of the calculations are compared with the existing experimental data. The effects of 
the density tails of the projectile nuclei as well as the role of shell admixtures and 
^ ' short-range correlations are analyzed. 

o ■ 

2 ■ 1 Introduction 

^ \ 

2 ' The study of the matter distributions in the Borromean nucleus ^He and the neutron rich 

isotope ^He is an actual problem discussed in many papers (see, e.g., in and refs. therein). 
Investigations give a reason to expect long tails (halo) of the neutron densities. The first 
result on the enhancement of the strong interaction radii Rj for ^'^He was deduced in j2j 
^ I from measurements of their interaction cross sections with different target nuclei at fixed 

^ ' energy 800 MeV/nucleon. The same conclusion was made in p] about the ^'®He rms-radii. 

Later in the model densities in the form of a sum of two Gaussian functions |5j were 
^ . fitted to the data and a conclusion for a significant neutron halo in the outer region of these 

nuclei was made. As to the differential elastic scattering cross sections, they were studied in 
collisions of ^'^He with protons, and the matter distributions of these nuclei were tested using 
different phenomenological and theoretical models. Elastic scattering angular distributions 
of ^He + p, ^He + ^He {Eiab = 151 MeV) and ^He + p, ^He + p {E < lOOA MeV) were 
calculated in [0] and [Zj, respectively, by means of the real part of the optical potentials 
obtained microscopically by using the realistic M3Y-Paris effective interaction jHl El and 
taking into account the neutron and proton densities of helium isotopes from |3] and also 
derived from the cluster-orbital shell- model approximation (COSMA) [5] , TT' , [T^ , . It was 
shown in |7j that elastic scattering is a good tool to distinguish between different density 
distributions. The results of [7j were compared also with those of the alpha-core approach 
with the fully nonlocal effective interactions ^Ij and the no-core model based on the large 
space shell- model calculations (LSSM) [IS] , [15] , [H] . Elastic scattering was recently studied 
also in jTB]. 

At the same time, small experimental information on the energy dependence of the total 
reaction cross sections of ^'^He with nuclei is available. Generally, from such data for the well- 
known densities of the target-nuclei, one can get additional constraints on the parameters of 
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the tested models of the halo nuclei. Below, we present our calculations of the total reaction 
cross-sections and compare them with the respective data on interactions of ^'^He with '^^Si 
at energies E/nucleon ~ 10^ 50 MeV 



2 Some formulae 



Following |2Tj we start with the expression for the total reaction cross-sections in the frame- 
work of the optical-limit model of the high-energy multiple-scattering theory 



oo 

aR = 271 J dbb (l-e-^(^)) . (2.1) 







Here, the eikonal phase for the scattering of the projectile nucleus (p) by the target nucleus 
{t) is given in by the expression 

x{b)=aNN j d^Spd'^St plisp) pl{st) fid, ^ = b - Sp + St. (2.2) 

For the isospin averaged total nucleon-nucleon cross-section a^N we use the parametrization 
from [22] • The vectors s, ^ are in the impact parameter b plane perpendicular to the z- 
axis along the relative momentum kj of the two scattered nuclei. The thickness densities are 
given as 

oo 

p%s) = I dz p°{^^^T^), (2.3) 



where the bare density distributions p°(r) of point nucleons in respective nuclei depend on 
their coordinates in the center-of-mass frames. The function /(^) is defined by the form of 
the nucleon-nucleon interaction amplitude. 

To present the phase integral in one-dimensional form, one uses the two-dimensional 
Fourier-Bessel transformation of all functions in the integrand 



oo 



where 



This results in: 



u{s) = ^ j e-''^^u{q)d\ = i- 1 Uqs)u{q)qdq, (2.4) 



oo 

u{q) = j e^^^ u{s)d^ s = 271 j Jo{qs) u{s)sds. (2.5) 







oo 



X{h) = ^NN^ j qdq Mqb) p%{q) p°t{q) /(g). (2.6) 



The original high-energy scattering theory uses only the transversal component q± = 
gsin'i9/2 of the total momentum transfer q = 2ksmd/2 and neglects its longitudinal com- 
ponent q\\ = qsm'd/2. Here q^ = q']_ + q^^-, and k is the momentum of the relative motion. 
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So, the theory is apphed to the scattering angles •& < ^J^lJkK^ where K is approximately the 
sum of the radii of both the nuclei^. So, in eq. ()2.6|) the use of g_L = g allows one to take the 
3-dimensional Fourier transforms of the functions p°(r) and /(^) instead of the 2-dimensional 
transforms of their profile functions and /(s). We remind that the form factors of the 

point p°(g) and nuclear p(g) densities obey the relation 

p(g)=P^(g)p°(g), (2.7) 

where the nucleon form factor pAr(g) is suggested to be the proton one and has a dipole form 
which is transformed to the Gaussian function at small q: 

pN{q) = (l + ^^Y" ^ exp(-gVo%_/6). (2.8) 

At the same time the form /(g) of the NN-scattering amplitude is known to be 

j{q) = exp(-gV^,_/6). (2.9) 

The squared values of the rms radii of the nucleon and the NN-interaction, known from 
the experimental data as Tq^^^ = 0.658 /m^ and rj^rms — 0-66 ± 0.03 /m^, almost coincide. 
It was concluded in [21] that a small difference between these rms radii does not influence 
visibly the calculated total reaction cross sections, and, thus, one can use approximately 
PN^q) = f{q)- Then taking into account Eqs. (|2.7|) - (j2.9p we get the phase (j2.6p in a simple 
form 

oo 

Xib) = a^N ^ j qdq Jo{qh) p°p{q) pt{q). (2.10) 



The form factors pt{q) can be calculated by using the data on nuclear density distributions. 
Here, we use the realistic density in the form of the symmetrized Fermi density distribution 

^^^^ ^ sinh R/ a 1 

cosh i?/a + cosh r/a 1 + exp[(r — i?)/a] ' 



and the respective form factor obtained in [2Z] (see also J8^) in the form 

singi? 



^ Ana d 

p{q) = -po — T 

q dq 



sinh naq 



3A 



(2.12) 



Here R and a are the radius and diffuseness parameters, respectively. Now, taking different 
models for the projectile nuclei we can test them by comparison of the calculated cross 
sections with the existing experimental data. Note that if for the densities of the projectile 
and the target nucleus one takes the Gaussian functions and uses (|2.9|) for /(g), then the 
integration in (|2.6|) can be performed explicitly [21], and one obtains the result of [221 
frequently used by experimentalists analyzing experimental data 



^ For estimations of the validity of this approximation one can see |26| . 
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where aj^ = (2/3)r^^^^. In the case of the zero-range approximation = 0. For the 
repulsive Coulomb field the trajectory of the incident nucleus defiects from the scattering 
center. This effect can be taken into account replacing he impact parameter b in the phase 
x{b) t by the distance of the closest approach be in the Coulomb field 

b ^ bc = a + Va^ + b\ (2.14) 

where a = ZpZte^ /2Ec,m. is the half-distance of the closest approach at b = 0. 

In our work, we use models for the projectile nuclei which are known in a table form and, 
thus, we apply expression (|2.1(jp in the calculations where the density Ptiq) of the target 
nucleus '^^Si is taken in a realistic form 1)2.111) with R = 3.085 fm and a = 0.563 fm from 
jSO|. The Coulomb distortion effect is also taken into account by using ()2.14|) . 

3 Density distributions and calculated cross section 

We apply the model for the s-p and s-d shell nuclei, developed in Ref. (M- model), to 
calculate the ^'^He densities. The densities of these projectile nuclei are used to study the 
effect of the Jastrow-type short range NN-correlations (SRC), on the total reaction cross- 
section ^'^He +^'^Si measured at different energies (THlEni- 

In the M-model, a general expression for the one-body density matrix of N = Z, s-p and 
s-d shell nuclei jHIj has been obtained using the factor cluster expansion method (see Ref. jH^I 
and Refs. therein). That expression depends on the harmonic oscillator (HO) parameter B 
(B = (^/(mcu))^/^), the occupation numbers Ani of the various states and the correlation 
parameter f3 that comes from the Jastrow type correlation function 

f{rij) = 1 - exp[-/5(ri - rj^], r^j = 1^ - rj] (3.15) 

which introduces short range correlations. The effect of correlations introduced by the func- 
tion /(r) becomes large when the correlation parameter f3 becomes small and vice versa. 

The expression of p(r) (in the two body approximation for the cluster expansion) has the 
form 

PW = A^[P5D(rB) + 02(rB,/5)], rB=r/B (3.16) 

where PsDi''^B), the uncorrelated density distribution corresponding to the Slater determinant 
HO total wave function, comes from the one-body term of the cluster expansion of the 
one-body density matrix and the term 02{rB,P) comes from the two-body term. Their 
expressions as well as the expression of the normalization factor N are given in 01^ . Similar 
expressions for the form factor and the momentum distribution have also been given in Ref. 

m 

The parameters B, j3 and Ani can be determined, for each nucleus, by fit to the experi- 
mental form factors and momentum distributions. As such data are not available for ^'^He 
we will use the data of the total reaction cross-section and the rms-radius. That is, we 
involve the parameters S, (3 and the four occupation number Ais, Aip, Aid and A2s in the 
fitting procedure. 

We used as experimental data the well known results of Tanihata et al |3] (T-model) 
for the density distributions obtained by fit the total reaction cross-section of ^'^He +^^C 
at 800 MeV in calculation with the model discussed in Sec. 1, using for the target nucleus 
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the Gaussian form Pg{i^) fitted to the reahstic density in the outer region and using for the 
projectile nucleus the harmonic-oscillator form 



Pxir) 



X -2 



.^g(r/fe)2 



(3.17) 



where X = Z,N , and a, h are the fitting parameters within the T-model. 

Expression ()3.17p and the corresponding expression of the rms-radius have been used as 
the "experimental" data of the projectile nuclei ®'^He. The values of the parameters B, (3 and 
Ani which have been found by a least squares fit of the theoretical neutron and proton density 
distributions (expression ()3.1fij) ) to the "experimental" ones with the constrain that the two 
densities have the same rms-radius are given in Table I. In the case of proton densities there 
is one free parameter, the HO parameter. For that reason the proton density is taken to be 
the same for both models. 



10 



10 



10 



10 































^ \ 


(a) 


\ " \ ^ 


\ \ 




. \ . \ 


3 2 


4 6 8 




r (fm) 



2500 



2000 



1500 



1000 



500 




10 20 30 40 

E/A (MeV) 



50 60 



Figure 1: (a) Point density distributions of He-6 calculated using the M-model [Ej. Dashed 
line: neutron density; dash-dotted line: proton density; bold curve is a bundle of eight curves 
(see the explanations in the text); (b) Total reaction cross section for He-6 + Si-28. The 
curves are as in a). 




Figure 2: The same as in Fig. 1 from calculations within the M-model J7] where only the 
parameter beta was used in the fitting procedure. 

In order to study the effect of SRC on the density and the total cross section we examined 
two cases. In the first case we have taken the neutron occupation numbers to be the ones 
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Table 1: The values of the parameters B, (3 and Ani , for the neutron and proton densities 
of ^'^He, obtained by fit of the expression ()3.16p to expression ()3.17|) with the constrain that 
the two densities to give the same rms-radius. 



Nucleus 


B[im\ 


/?[fm]-2 


Ais 


Alp 


Aid 


A2s 




6He(n) 


1.8249 


6.0001 


1.9999 


1.0889 


0. 


0.91124 


2.5865 


8He(n) 


1.8075 


5.1558 


2.0000 


2.9056 


0. 


1.0944 


2.6587 


''He(p) 


1.53 


oo 


2. 


0. 


0. 


0. 


1.72 


8He(p) 


1.53 


oo 


2. 


0. 


0. 


0. 


1.76 



given in Table I, we have given to the correlation parameter (3 various values from /5 = 4 fm~^ 
to = 50 fm~^ and for each value of (3 we determined the parameter B so that the neutron 
rms-radius to be the same with the one obtained within the T-model. In figure 1 we present 
the point density distribution for ^He calculated in this case. The bold curve is the bundle 
of eight figures. It is seen that the neutron densities of the M-model have lower tails at 
r > 5fm compared to the one of the T-model. At the same time, we see that the cross- 
sections are in good coincidence. This means that the main contribution comes from the 
region r = ^ 5fm, i.e., up to r equal to about two rms radius of the ^He nucleus (the 
rms-radius of the mater distribution is 2.331 fm for both model). From a physical point of 
view, the increase of the SRC at small (3 enhances the nuclear radius, but in the considered 
case one can think that this effect is compensated by reducing the HO parameter B which 
results in "squeezing" the wave functions of the basis. For that reason we considered the 
second case where the calculations were made for various values of the parameter (3 in the 
range 4 to 50fm~^ while the HO parameter B had the value 1.8249 fm and the parameters 
Ani were taken as before from Table I. The results are presented in Fig. 2. The calculations 
in Ref. PI] showed that the changes of the values of (3 in that region give very significant 
changes in the form factors and the momentum distributions of nuclei. In the present paper, 
we study the total cross-sections and as a result, we obtained the corresponding spread of 
the rms radius of the neutron density for ^He in the limit 2.6030 ^ 2.5757 fm. However, this 
does not noticeably influence the cross-sections. Once again, the cross-sections of the M- 
and the T- models are in coincidence, while the rms radii are slightly different. 

The same conclusions can be made for ^He nucleus from Figs. 3 and 4, where the 
results of the calculations within the M-model are presented. In Fig. 3, one can see the 
bunch of nine curves for the neutron densities when the pairs of the parameters (3 and B 
which give the rms-radius Rrms,n = 2.6587 fm are in the regions /5 = 3.1588 -r- 50fm~^ and 
B = 1.7850 -i- 1.8261 fm. The occupation numbers are the ones given in Table I. For all 
the neutron densities the rms-radii are in coincidence with the T-model Rrms,n = 2.6857 fm. 
The matter density radius is also the same in both models Rrms,A = 2.6857 fm. The values 
of the parameters of the curves in Fig. 4 are: f3 = 3.1558 — 10.1558 fm~^, B = 1.8075 fm and 
Ani are the same as in Fig. 3. The obtained values of the rms radii of the neutron densities 
are in the region 2.7172 -h 2.6676 fm. We see that in both figures the slopes of the densities 
from the M-model begin to differ from that of the T-model at about 5 fm, where the densities 
fall down by about two orders of magnitude comparing to the values in the nuclear center. 
Their difference in the range of twice the rms-radius does not affect the total cross-sections. 
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Figure 3: (a) Point density distributions of He-8 calculated using the M-model [17]. Dashed 
line: neutron density; dash-dotted line: proton density; bold curve is a bundle of neutron 
densities; (b) Total reaction cross section for He-8+Si-28. 



4 Discussion and conclusions 

One should mention that in the M-model there were no calculations of the binding ener- 
gies and the structure of excited states of these nuclei. However, large-scale shell-model 
calculations (LSSM) of the wave functions and characteristics of light exotic nuclei, such 
as ^'^He and ^'^^Li, were reported by several groups, e.g. [TTl ESI IMl IMl EEl- They are 
usually performed by using interactions obtained directly from NN g-matrices with various 
NN-interactions as their base. In E3S|, the wave functions for ^'^He were calculated within a 
complete 4hw model space by using the g- matrix with interaction from . All calculations 
were performed by the shell-model code OXBASH [38j. For ^ife the harmonic-oscillator 
single-particle wave functions were used, while for ^'^He they were the Woods-Saxon ones. 
It was concluded in |36j that ^He corresponds to halo nucleus, while ^He is a neutron skin 
nucleus. In [HE] (and refs. therein), the K-harmonic method was developed and calcula- 
tions were based on the assumption of the lp3/2 state for each of the valence neutrons related 
to the alpha-core center (COSMA- model). The admixtures of the 3-body forces is included 
to get the proper value of the binding energy. 




Figure 4: The same as in Fig.3 with differences in the fitting procedure (see the explanations 
in the text). 
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In Fig. 5, we show the matter density distributions for ^He obtained from the calcula- 
tions in the framework of the COSMA-model ^3] (dash-dots) and the LSSM-model 
(dots). The M- and T-model densities are shown by the solid bunch and the dashed curve, 
respectively. First of all, one can see that in the outer region the COSMA- and LSSM-model 
densities have extended tails which are followed by exponential asymptotics. Here, the dot- 
ted curve exceeds the dash-dotted one which results in the enhancement of the corresponding 
cross section. It is seen that both the cross sections are about 200 500 mb above those 
calculated by the M- and T-models and the experimental data. The respective rms matter 
radii are Rrms,N = 2.560 fm (for the dash-dotted curve) and Rrms,N = 2.956 fm (for the 
dotted curve). In this connection we refer to [SH] where the coupling of elastic and elastic 
breakup channels were shown to play an important role in processes with weakly bound 
nuclei, and this can diminish the cross section obtained without accounting for the addi- 
tional channels. Recently, the theory of the coupled channels was developed in ^01 to 
study the breakup reactions of halo nuclei. In general, the implied optical-limit model needs 
higher order NN-scattering terms to be included as had been done in the multiple-scattering 
theory [^I^. For example, it was shown in ^Sl that these terms provide a negative inter- 
ference and thus diminish the p + ^He differential elastic cross section calculated by taking 
into account only the single-scattering term. The in-medium effect, which changes the free 
NN-scattering cross section aNN in the phase ()2.2p can be included as well. As shown qual- 
itatively in [21], this effect can decrease the total nucleus-nucleus reaction cross sections by 
about 10%. All the mentioned problems regarding the mechanism of the nucleus-nucleus in- 
teraction are still under investigations of the current theory. However, our study of the total 
reaction cross-sections of the ^'^He interactions with the "^^Si target leads to some general 
conclusions which, we hope, will not be affected by a possible further modification of the 
theory. 




Figure 5: (a) Matter density distributions for He-6 calculated with the COSMA model 
[21,22] (dash-dotted hue) and the LSSM [27] (dotted hue). The M- and T-model densities 
are given by solid bunch and dashed lines, respectively, (b) Total reaction cross section for 
He-6+ Si-28. The curves are as in (a). 

The following conclusions can be drawn from our results: 

• The rms radii are more sensitive nuclear characteristics than the total reaction cross 
section with respect to changes of the density distributions. The radius depends 
strongly on variations of the density in the nuclear interior while the behaviour of the 
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total reaction cross-section depends mainly on the peculiarities of the nuclear surface 
and weaker on the central nuclear region where strong absorption takes place. 

• The previous item explains why variations of the densities due to SRC do not affect 
the total cross sections while they are significant in calculations of the form factors and 
of the single-particle momentum distributions in nuclei at large momentum transfers. 

• The fitting procedure in our work showed that there is not a Ici-shell contribution to 
the wave functions of ^He and ^He. 

• In the case of the harmonic-oscillator basis, the Gaussian type of the wave function 
asymptotics does not affect significantly the cross sections and the latter are formed 
by taking into account the interior region of the projectile nucleus up to about only 
two rms nuclear radii. Otherwise, if one gets a basis for realistic wave functions, which 
have an exponential asymptotics, the region of the main contribution becomes larger 
and achieves the value of about four and larger times the rms nuclear radius. 
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